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We analyze transport of local magnetization and develop schemes to control transport behavior 
in finite spin-1/2 Heisenberg chains and spin-1/2 Heisenberg two-leg ladders at zero temperature. 
By adjusting parameters in the Hamiltonians, these quantum systems may show both integrable 
and chaotic limits. We provide examples of chaotic systems leading to diffusive and to ballistic 
transport. In addition, methods of coherent quantum control to induce a transition from diffusive 
to ballistic transport are proposed. 
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I. INTRODUCTION 

A complete understanding of transport behavior in 
I many-body systems is one of the utmost challenges in 
fundamental studies of non-equilibrium statistical me- 
chanics. In the classical domain, it is widely believed that 
chaotic systems should show diffusive (normal) transport, 
whereas integrability should be associated with ballis- 
tic (abnormal) transport [Ij, although normal trans- 
port has also been verified for a non-chaotic system 
In the quantum domain, the conditions that determine 
specific transport behaviors are still under debate, but 
here also, the main conjecture is of the correspondence 
integrable-ballistic and chaotic-diffusive Several the- 
oretical approaches have been undertaken to address this 
issue, including: attempts to derive the Fourier law from 
a microscopic foundation by applying the Hilbert space 
\ average method Q and by numerically studying trans- 
port of heat in finite chaotic and non-chaotic systems 
coupled to heat reservoirs [1, Q; analysis of the trans- 
port behavior of local magnetization in isolated finite 
systems at zero temperature [sl; new advances toward 
the problem of quantum thermalization Q ; comparisons 
of the results for conductivity in integrable and chaotic 
systems at finite tem per ature and in the thermodynamic 
limit [13, El, [12, a subject of intense discussion 

here being the possibility of ballistic transport in non- 
integrable quasi-one-dimensional systems [isl. [l^ [l7| : nu- 
merical studies of spin diffusion at long times at infinity 
temperature [Tsj : as well as studies of transport proper- 
ties near the metal-insulator transition [19, 20, 2li|. 
Investigations of transport behavior in the particu- 
' lar case of quasi-one-dimensional spin-1/2 systems have 
been highly motivated by experiments in low-dimensional 
magnetic compounds, such as copper oxide (cuprate) sys- 
tems, where ballistic behavior has been observed for heat 
conduction [S^, |2^, and also for magnetization, as 
revealed by nuclear magnetic resonance (NMR) experi- 
ments (25| . These compounds are described by models 
of interacting spins-1/2 arranged in structures such as 
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chains, two-leg ladders, and square lattices [26[. Clean 
spin-1/2 Heisenberg chains with only nearest-neighbor 
interactions are integrable models solved with the Bethe 
Ansatz method [l^], whereas two-dimensional lattices 
are chaotic [1^ and two-leg spin ladder systems become 
chaotic when the interchain and intrachain interactions 
are of the same order [2§| . 

In Ref. [1], an isolated isotropic finite spin-1/2 Heisen- 
berg chain at zero temperature and with only nearest- 
neighbor interactions was considered in the analysis of 
transport of local magnetization in both cases: when the 
system was clean and therefore integrable, and when ran- 
dom on-site disorder led to the onset of quantum chaos. 
Free boundary conditions were assumed. A bouncing be- 
havior of the local magnetization was observed in the 
integrable regime and interpreted as a hint of ballistic 
transport, whereas for the chaotic system the local mag- 
netization showed an exponential relaxation to equilib- 
rium, which was considered a reflection of diffusive trans- 
port. 

The first part of this paper is also dedicated to the 
investigation of transport of local magnetization in spin- 
1/2 Heisenberg systems, but different ways to induce 
chaos are dealt with. The goal is to resolve whether 
non-integrability may have a unique correspondence with 
the exponential decay verified in The integrability- 
breaking terms considered are: on-site disorder described 
either by (i) randomly distributed Zeeman energies (soj 
or by (ii) a few defects placed on specific positions of 
the chain [sij : and couplings to more surrounding spins 
such as (iii) next-nearest neighbor interactions and (iv) 
interchain interactions, as typical of two-leg spin ladder 
systems. An exponential decay is observed only for (i) 
and (iv). For the finite systems studied, few defects 
lead to a behavior that sometimes resembles localiza- 
tion, the exponential behavior appearing only under very 
special conditions, while next-nearest neighbor interac- 
tions generate oscillations of the local magnetization that 
are similar to those seen for the integrable system, al- 
though faster. These observations support the conjecture 
that chaos might not be a sufficient condition for normal 
transport. 

The second part of this work focuses on the analysis 
of methods of coherent quantum control, the so-called 
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dynamical decoupling (DD) schemes, as potential tools 
to manipulate transport behavior. DD schemes consist 
of sequences of external control operations that average 
out unwanted contributions to the system Hamiltonian. 
These methods have long been applied in NMR spec- 
troscopy [13, [s^l , where the goal is to modify the nuclear 
spin Hamiltonian to eliminate or scale selected internal 
interactions. More recently, DD has addressed the re- 
moval of interactions between the target system and its 
environment [3^ and has been put in a general control- 
theoretic framework [35|. It has also been considered in 
studies of transport of information [s^. Here, we in- 
troduce DD sequences that suppress the effects of terms 
leading to quantum chaos, the purpose being the ap- 
proach to the ballistic transport behavior verified for the 
integrable system. 

The paper is organized as follows. Sec. H explains how 
the identification of the chaotic regime is performed and 
describes the models to be considered. Sec. HI presents 
the results for the transport of local magnetization for 
systems in different regimes. Sec. IV introduces DD se- 
quences to control transport behavior and shows, as an 
illustration, results for the sequence that cancels the ef- 
fects of on-site disorder. Concluding remarks are given 
in Sec. V. 



II. QUANTUM CHAOS AND SYSTEM MODEL 

A. Signature of Quantum Chaos 

For quantum systems, chaos may be identified by ana- 
lyzing the distribution of spacings s between neighboring 
energy levels [13, HI] . Quantum levels of integrable sys- 
tems tend to cluster and are not prohibited from crossing, 
the typical distribution is Poissonian: 



Pp{s) = exp(-s). 



(1) 



In contrast, non- integrable systems show levels that are 
correlated and crossings are strongly resisted, the level 
statistics is given by the Wigner-Dyson distribution. The 
exact form of the distribution depends on the symmetry 
properties of the Hamiltonian. In the case of systems 
with time reversal invariance it is given by: 



Pwd{s) = 7rs/2exp(-7rs^/4). 



(2) 



To analyze the transition from integrability to chaos, 
the quantity 77, defined as 
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tem 77 ^ 1, while for a chaotic system 77 — >• 0. The criti- 
cal value below which the system is considered chaotic is 
chosen to be 77 = 0.3 [40l |. 



B. Heisenberg Model 

We study homogeneous and isotropic spin-1/2 Heisen- 
berg chains with open boundary conditions, as described 
by the Hamiltonian: 



H — Hz + Hnn + HjMNN 

L L-1 L-2 

= ^nS'n + JSn-Sn+1 + ^ J' 8^-8^+2 -(4) 

71—1 n—1 n—1 

Above, h is set equal to 1; Sn — (?n/2 is the spin op- 
erator at site 71, CT^^'^'^ being the Pauli operators; and 
L corresponds to the number of sites. The parameter 
LDn is the Zeeman splitting of spin n as determined by a 
static magnetic field in the z direction. The system is 
clean whenever all sites have the same energy splitting 
LOn = OJ, and it is disordered when defects characterized 
by different energy splittings ojn = ^ + dn are present. J 
and J' are the interaction strengths of nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) couplings, re- 
spectively, and are assumed to be constant. 

All calculations in this work are performed in the basis 
consisting of eigenstates of the total spin operator in the z 
direction, S"" = X^Li ^n- In this basis, the NN and NNN 
Ising interactions, S^S^^i and S^S^j^,2, contribute to the 
diagonal elements of the Hamiltonian, while the XY- 
terms, S!^^S:ll_^i + S'^S'^_,_]^ and S!^S^^^2 + consti- 
tute the off-diagonal elements. The role of the XF-terms 
is to transfer excitations through the system by exchang- 
ing the position of nearest and next-nearest neighboring 
spins pointing in opposite directions. 

Also considered here are two-leg spin ladder systems 
corresponding to two coupled spin chains as described 

by 



L/2 L/2-1 

^ ^ ^n,mSj^ „^ -f- ^ ^ J Sn,m-Sn+l,ri 
n—1 n—1 



E 

m— 1 
L/2 

J±Sn,l-S. 



(5) 



was introduced in [39|, where sq ~ 0.4729 is the first in- 
tersection point of Pp and Pw d ■ For an integrable sys- 



where ujn and J are as before, tt^ determines the chain 
in which the site is positioned, and J± characterizes the 
strength of the interchain interaction. 

In order to derive meaningful level spacing distribu- 
tions, before diagonalizing the Hamiltonian and unfold- 
ing the spectrum [stI . [38i] , all trivial symmetries of the 
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system need to be identified. The analysis of symme- 
tries is necessary, because a Poisson distribution may ap- 
pear whenever Wigner-Dyson distributions from different 
symmetry sectors are mixed, which may lead to erro- 
neously interpret the system as integrable. In both mod- 
els considered here, is conserved, therefore, instead 
of diagonalizing matrices of dimension 2^ , we study the 
largest subspace. For L even it corresponds to the sector 
with ^ and dimension N = (^^J = L\/[{L/2)\]^. 
Depending on the parameter values, Hamiltonians (|4l) 
and ^ may also exhibit the following symmetries [28|: 
invariance under lattice reflection, which leads to par- 
ity conservation; and conservation of total spin S'^ = 
that is [H,S^] = {S^ symmetry). Notice 
also that the Heisenberg model with a magnetic field does 
not commute with the conventional time-reversal opera- 
tor, however the distribution associated with its chaotic 
regime is still given by Eq. ([2]) , as discussed in [H, [s^l • 



III. TRANSPORT OF LOCAL 
MAGNETIZATION 



In studies of transport properties, the most popularly 
used method is the Green-Kubo formula [4l[. However, 
the application of this formula for the treatment of heat 
transport has been criticized and the use of the Hilbert 
space average method to demonstrate the emergence of 
heat diffusion from microscopic models has been sug- 
gested as an alternative 0. This approach has been 
extended to the analysis of transport of magnetization 
in Ref. [S]. 

Here, as in Q , we study the transport of local magne- 
tization as defined by 



L/2 



(6) 



n=l 



where \ip{t)) is the state of the system at instant t written 
in the basis of . The initial states considered come from 
the sector S*^ = 0. For system (U)), \ip{0)) has all spins 
pointing up placed in the first half of the chain, while 
the remaining down-spins appear in the other half. For 
the two-leg system, the initial state has all up-spins in 
one chain and all down-spins in the other. We assume 
L = 12, which leads to Af(0) = 3 in both cases. 

Integrable system. The clean system with only nearest- 
neighbor interaction, as described by H ^ with dn — 
and J' = 0, corresponds to an integrable model solved 
with the Bethe Ansatz method 27]. The dynamics for 
the local magnetization is shown on the left panel at the 
top of Fig. [TJ The bouncing behavior suggests ballistic 
transport [8]. In what follows, we compare this result 
with the time evolution of local magnetization for chaotic 
regimes induced by different ways. 



A. Chaos induced by on-site disorder 

In the chain with only nearest-neighbor interactions, 
as given by H ^ with J' = 0, chaos may be induced if 
one or more defects are present. 

Random on-site disorder. Assume that the Zeeman 
energies are given by ujn = uj + dn, where dn are uncor- 
related random numbers with a Gaussian distribution: 
(dn) = and {dnd„i) = d'^Sn,m [30| . The transition from 
integrability (d = 0, rj ^ 1), to chaos (0.05< d/J < 0.7, 
rj < 0.3) is indicated by 77, which is computed in the sec- 
tor 5'^ = and is shown on the right panel at the top of 
Fig. [1] Notice that as d/J — > 0, conservation of parity 
and total spin start playing a role. At d/J = 0, the cor- 
rect evaluation of 77 would need to take these symmetries 
into account 1281. 









' 1 ■ 1 








I.I.I 





200 


400 600 



Jt 



FIG. 1: (Color online.) Transport of local magnetization in 
a Heisenberg cliain witli only nearest- neighbor interactions 
described by © with J' = and L = 12. The value of 
UJ is irrelevant for the dynamics. Left top panel: M{t) for 
d„ — 0. Right top panel: rj computed in the sector — 
for Gaussian random on-site disorder, {d„) — and (dndm) = 
d^5„,„. Left bottom panel: M(f) for d/J = 0.05, 0.1, 0.15, 0.2, 
and 0.25 from bottom to top. Average over 20 realizations. 
Right bottom panel: M{t) for a sample of realizations with 
d/J = 0.15. 

On the left panel at the bottom of Fig. [U we show the 
time evolution of M{t) averaged over 20 realizations for 
five different values of d/J in the chaotic region. As we 
approach chaoticity, for d/J = 0.05 and 0.1, oscillations 
are still seen; whereas for d/ J = 0.15, 0.2, and 0.3, an ex- 
ponential decay of M{t) takes place reaching final values 
between 1 and 2. After the decay, the larger probability 
to find spins pointing up in the first half of the chain 
is reflected by the positive values of AI{t), which are ob- 
tained with the majority of the realizations. On the right 
panel at the bottom, we show the behavior of a sample 
of realizations with d/J = 0.15: for some of them M 
decays to equilibrium, M ~ 0, indicating an equal prob- 
ability to find up-spins in both halves of the chain, but 
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for the majority of the realizations M{t) remains positive 
throughout; hardly any curve reaches negative values of 
M. 

One defect. A single defect in the middle of a chain, 
at n = L/2 (or equivalently at n = L/2 + 1), may lead to 
quantum chaos. For L = 12, this happens when 0.15 < 
d-e/J ^ 2.0, as discussed in [31]. The transition to chaos 
is shown by rj on the left panel at the top of Fig. [51 The 
curve is obtained in the 5^ = sector, but, as mentioned 
before, for d^/J ~ the correct evaluation of rj would 
require also the consideration of the S'^ symmetry and 
parity conservation j2^ . A way to break the symmetry 
and deal with larger subspaces even when the system is 
integrable consists of including defects at the edges of the 
chain [3l|,|4|. 
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FIG. 2: (Color online.) Transport of local magnetization in 
a Heisenberg chain with only nearest-neighbor interactions 
described by // Q with J' = and L = 12. Left panels: one 
defect on site 6, d„ =0 for n 7^ 6. Right panels: two equal 
defects on sites 6 and 7, d„ = for n 6,7. Top panels: 
dependence of r; on the value of the defect(s). Bottom panels: 
M{t) for chaotic systems. 

The transport of local magnetization is shown on the 
left panel at the bottom. For de/J — 0.3, 0.5, and 1.0, 
even though rj indicates chaoticity, an exponential de- 
cay is not observed. Instead, for the two smaller values, 
especially for the smallest one, partial revivals are ver- 
ified, whereas for d^/J =1.0 there occurs localization 
of the up-spins in the first half of the chain. In con- 
trast, a defect placed on site 7, which shows exactly the 
same behavior for 77, does lead to an exponential decay 
of M(t) - see right panel of Fig. [3] [43]. The reason for 
the different results is the following. The initial state 
I T1T2T3T4 jVlVtr i8i9iioiiii-i2) is directly coupled only 

with the state | Tit2T3T4 hleh i8i9iioJ.iiii2>- When 

the defect is on site 6, the hopping of the up-spin from 
site 6 to site 7 is not favorable, since the state loses the 
extra energy from the defect site and the positive Ising 
energy coming from the pair of parallel spins tsTe- This 



explains why it becomes easy to localize the initial state 
by increasing the defect value. Contrary to that, if the 
defect is on site 7, the Ising energy lost by breaking the 
pair of parallel spins is regained by placing the up-spin 
on the defect. In this scenario, directly coupled states 
may be very close in energy, which favors delocalization. 

Two defects. Two equal defects in the middle of the 
chain, at n = L/2,L/2 -f 1, may also lead to quantum 
chaos. For L = 12 this happens when 0.15 < de^r/J < 
2.0. The transition to chaos is shown by 77 on the right 
panel at the top panel of Fig. O The curve is obtained 
by taking both symmetries into account: conservation of 
parity and Sz ■ The relaxation of the local magnetization 
for d^j/ J — 0.3, 0.5 and 0.65 is shown on the right panel 
at the bottom. The decay is not exponential and M = 
is never reached, most up-spins tending to localize in the 
first half of the chain. This may again be understood 
by comparing the energies of the initial state and of the 
state it is directly coupled to. 

The whole spectrum is required to obtain the plots for 
■q presented in Figs. [T] and [21 therefore the decision to 
deal with relatively small systems, L = 12. This choice 
was a good compromise leading to sufficient statistics and 
the possibility to run various realizations of random on- 
site energies. In addition, notice that the behaviors of 
the local magnetization obtained with L = 12 are also 
reproduced with 10 and 14 spins, as shown in the two 
panels of Fig. [3l Thus, for the analysis developed in this 
paper, a system with 12 spins is sufficiently adequate. 
In order to simulate the time evolution of M{t) in much 
larger systems, we could resort, for example, to the very 
efficient algorithm recently proposed in Ref. [i^ . 
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FIG. 3: (Color online.) Transport of local magnetization in 
a Heisenberg chain of different sizes described by _ff Q with 
J' = 0. Left panels: clean chains, dn = for 1 < n < L. 
Right panels: chaotic systems with one defect in site L/2-I- 1, 
d„ = forn / L/2 + 1 and dL/2+ilJ = 0.5. 



B. Chaos induced by additional interactions 

In a clean Heisenberg system with dn — 0, chaos may 
be induced by adding frustrating next-nearest neighbor 
interactions where J' ~ J or by keeping J' = and 
adding interchain interactions where J± ^ J [1^. In 
both cases, total spin, total spin in the z direction, and 
parity are conserved. 
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Next-nearest neighbor interactions. Hamiltonian (|4]) 
with ci„ = and J = J' describes a chaotic system j29i] . 
However, the transport of local magnetization shows os- 
cillations similar to those observed for the integrable sys- 
tem, although at a faster rate, as seen on the left panel 
of Fig.m Therefore, if the bouncing behavior of M{t) in 
isolated systems at zero temperature is indeed a signa- 
ture of ballistic transport, integrability is not a necessary 
condition for abnormal conductivity. 

Interchain interaction. Hamiltonian ([5]) with dn = 
and J — J± also describes a chaotic system [2^. In this 
case, as shown on the right panel of Fig. 01 a very fast ex- 
ponential decay of the local magnetization to equilibrium 
is indeed verified. 
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FIG. 4: (Color online.) Transport of local magnetization in 
a clean Heisenberg system. Left panel: H ((4]) with L — 12, 
d„ = 0, and J' = J. Right panel: H © with L = 12, d„ = 0, 
and J± = J. 



The distinct behaviors of the transport of local mag- 
netization verified for different chaotic systems ~ Figs.[Tl 
[21 [HI and [D ~ prevent a clear correspondence between 
chaoticity and diffusive transport. New light may be shed 
to the problem if this analysis is extended to the trans- 
port of other quantities, such as heat in open systems. 



rection perpendicular to z, for example X ^ clS determined 
by the operator. 



Px = exp -iTT^ I = exp(-i7rS"^). 



The interaction terms remain undisturbed, but the se- 
quence of these rotations leads to the cancellation of 
the one-body terms at every t2p — pTc, p 6 N, where 
Tc = 2At is the cycle time (for details about rele- 
vant frames and pulse generation see [i^ and references 
therein). The propagator at a time t2p is then given by 



U{pT,) ^ P^U{t2p,t2p-i)Px . ■ ■U{t2M)PxU{ti,Q), (7) 

where U{t) — Texp[— i j^Hdu] and T denotes time or- 
dering. By adopting the notation 

U+ = exp[-i(i7;, + HNN)^t] = exp[-ii7iAi], 



U- = PxU{tj,tj^i)Px= Px{PxPl)U{tj,t,^i)Px 



-lexp 



+ HNN)e~"'^'' ] At 



= — cxp[~i{~Hz + HNN)At] — — exp[~iH2At], 

where 1 is the identity operator, we rewrite the propaga- 
tor as 

U{pTc) ^U-U+...U-U+ ^ exp[^iHpTc]. 

Above^ H = E^o^^''-' is the average Hamiltonian 
[s^ . [33| , and the terms in the sum are obtained by using 
the Baker-Campbell-Hausdorff expansion. The first two 
terms are 



IV. CONTROL OF TRANSPORT BEHAVIOR 

We propose to control transport behavior by applying 
DD methods. DD schemes consist of sequences of exter- 
nal control operations that average out unwanted terms 
of the system Hamiltonian. In the case of spin systems, 
the control operations correspond to very strong mag- 
netic fields (pulses) able to rotate the spins and time- 
reverse the system evolution [H, US]- Here, we assume 
the ideal scenario, where the pulses are arbitrarily strong 
and capable of performing instantaneous rotations, the 
so-called bang-bang controls Our goal is to elimi- 

nate the effects of the terms leading to quantum chaos 
and approach the transport behavior verified for the in- 
tegrable system shown Fig. [1] This can be achieved with 
the sequences described below, where each sequence han- 
dles a particular integrability-breaking term. 

On-site disorder. The effects of on-site disorder may 
be eliminated by rotating all spins after every interval of 
free evolution t^+i —tj = At, j G N, by 180° around a di- 



= —iHl+H2) = HNN, 

The sequence of pulses reshapes the Hamiltonian. In 
the ideal limit of At — > one recovers the Hamiltonian 
for a clean Heisenberg model with only nearest-neighbor 
couphngs, H sa Hnn, as desired. In Fig. [5l we show 
the transport of local magnetization as modified by the 
above sequence for two cases of on-site disorder: the 
top panels correspond to random defects and the bot- 
tom panels to two defects in the middle of the chain. 
The transport behavior for both situations - integrable 
system and disordered Heisenberg chain subjected to the 
DD sequence - closely coincide when the intervals be- 
tween the pulses are smaller than the reciprocal inter- 
action strength. At < J~^, as seen on the left panels. 
For short time evolutions, good agreement between the 
curves is still verified for At ~ J~^, whereas at longer 



6 




FIG. 5: (Color online.) Transport of local magnetization in a 
Heisenberg chain with only nearest-neighbor interactions de- 
scribed hy H ^ with J' = and L = 12. (Blue) Curves 
showing a fast decay to M{t) ~ 1 correspond to disordered 
systems in the absence of pulses. (Green) Bouncing curves 
represent disordered systems subjected to the DD sequence 
(O; they closely coincide with the (red) bouncing curves ob- 
tained with the clean Heisenberg system. Top panels: Gaus- 
sian random on-site disorder, {d„) — and {d„dm) = d?'5n,m 
with d/J = 0.2. Average over 20 realizations. Bottom pan- 
els: d(i,7 / J = 0.65 and d„ = for n ^ 6,7. Left panels: 
At/ J = 0.25. Right panels: At/ J = 1.0. Data is acquired 
after every Tc = 2At. 



times the accumulation of residual averaging errors be- 
come significant, this being more perceptible in the bot- 
tom right panel. To slow down error accumulation, ran- 
domized schemes as developed in [45,, 46] may be incor- 
porated to the pulse sequence. 

Interchain interactions. To decouple the two interact- 
ing chains in a two-leg spin ladder system, we may apply 
a sequence of 7r-pulses that rotates all spins of just one 
of the chains. The pulses do not affect the intrachain 
interactions, but frequently change the sign of the inter- 
chain interactions. To eliminate the interchain couplings 
in the three directions, we need a sequence of at least 
four pulses with two alternating directions, such as 

UiT, = 4At) = Py^2U{/At)P^^iU{M)Py^2U{M)P^sU{M), 
where 



L/2 



=exp -iTT^S-i^i , 



L/2 



p, 



2 = exp -ITT 



V 

n,2 



The chain which is subjected to the pulses is also al- 
ternated to guarantee the cancellation of the one-body 



terms as well. At any time t4p = pTc = ApAt, we then 
have H = Hnn,i + Hnn,2 + 0{At). 

Next-nearest neighbor interactions. A possible se- 
quence to eliminate next-nearest neighbor interactions 
involves eight 7r-pulses. The controls are applied to 
specifically chosen spins and are only viable if means ex- 
ist to distinguish them. The suggested pulse sequence 



U{Tc = 8At) ^ PsU{At) . . . P2U{At)PiUiAt), 
where 



Pi 



Pi 



P^ 



Pa 



P7 



Pr 



L(i-i)/4J 

n ^ 

k=0 
L(L-3)/4j 



L(L-2)/4J 

+ 4fc g~i''^^2 + ik ^ 

L(L-4)/4J 



fc=0 
L(L-2)/4j 



fe=0 
L(L-3)/4J 



fc=0 

L(i-i)/4J 



fe=0 
L(L-4)/4J 



fc=0 



A:=0 



At any time = pTc ~ 8pAt, we obtain H = Hnn/2 + 
0{At). Apart from a rescaling factor 1/2, we recover 
the Hamiltonian for a clean Heisenberg model with only 
nearest neighbor interactions up to first order in At. No- 
tice that the sequence cancels both next-nearest neighbor 
interactions and one-body terms. The Hamiltonians for 
the intervals of free evolution are given in the appendix. 
Since they no longer conserve Sz, simulations involving 
this pulse sequence need to consider the whole Hilbert 
space of dimension 2^. 



V. CONCLUSION 

The purpose of this work was twofold: to contribute 
to the ongoing discussion about transport properties of 
quantum many-body systems and to study the possibility 
of controlling transport behavior by resorting to dynam- 
ical decoupling methods. 

Our conclusions are based on the analysis of the trans- 
port of local magnetization in isolated finite spin- 1/2 sys- 
tems with free boundary conditions and at zero tempera- 
ture. Diffusive transport is associated with the exponen- 
tial decay of the local magnetization and ballistic trans- 
port with its bouncing behavior. Under these conditions, 
it was shown that a one to one correspondence between 
quantum chaos and normal transport does not necessar- 
ily hold; instead, different integrability-breaking terms 
may lead to both ballistic and diffusive behavior. 

Dynamical decoupling sequences capable of suppress- 
ing the terms of the Hamiltonian leading to quantum 
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chaos were proposed. The goal was to approach the bal- 
hstic transport behavior of local magnetization obtained 
with the integrable system. This was achieved by apply- 
ing pulses separated by intervals smaller than the recip- 
rocal interaction strength, as illustrated for two cases of 
on-site disorder: random defects and two defects in the 
middle of the chain. 

The manipulation of transport behavior in systems 
subjected to dynamical decoupling methods is a topic 
that deserves further exploitation. In particular, the 
analysis of heat transport in systems coupled to two heat 
reservoirs at different temperature is one of our future 
goals. Also important is the identification of real sys- 
tems where these ideas could be tested experimentally, a 



possible good candidate being a crystal of fluorapatite, 
as proposed in studies of information transport [3a | . 
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APPENDIX A: SUPPRESSING NEXT-NEAREST 
NEIGHBOR INTERACTIONS 

The pulses from Eq. dH) lead to the propagator 



U(T^) = PsU{At)PrU{At)PGU{At)Pr^U{At)PiU{At)P3U{At)P2U{At)PiU{At) 
= jPsP?. ■ ■ PiHP? ■ ■ ■ Pi)^U{At){Pj ...Pi) .. . Plu{At)Pi U{At\ 

= 1 exp{~iHsAt) . . . exp{-iH2At) exp{-iHiAt), 

I 

where, for L even and by using the notation, L(i-_3)/4J L(i-4)/4j 

ZC= 2^ W3+4fcS'|+4j, zd= 2^ UJ4+4kSl_^_4k 
k=0 k=0 

L(L-1)/4J L(i-2)/4J 



za 



k=0 






k=a 




1 








Xlo 


L/2 
k=l 


nx 
-1^2k 


Xle 


L/2-1 

\ ^ QX QX 

— 2-^ '-'2fe'-'2fe+l 
k=\ 


X2 = 


L-2 

\ ^ qx qx 

- ^n^n+2 
n=l 




Ylo 


L/2 
k=l 


qy 
^l'^2k 


Yle ■ 


L/2-1 

_ qy qy 
~ '^2k'^2k+l 
k=l 


Y2 = 


L-2 

\^ qy qy 

/ , '^n'^n+2 
n=l 




Zlo 


L/2 
k=l 


-l'-'2fe 


Zle 


L/2-1 

\ ^ qz qz 

- 2^ ^2k'~'2k+l 
k=l 


Z2 = 


L-2 

\ ^ qz qz 

/ , ^n^n+2^ 

n=l 



the Hamiltonians during the intervals of free evolutions 
are 



Hi = 


+za + 


zb + 


zc + 


zd H 


-XI0 + 


Xle4 


-YI0 + 


yie4 


ZI0 + 


Zle 


+ X2 


+ Y2 + 


Z2 


H2 = 


—za — 


zb + 


zc + 


zd H 


-XI0 + 


Xle 4 


Flo- 


yie4 


Zlo- 


Zle 


+ X2 


-Y2- 


Z2 


H3 = 


—za — 


zb- 


zc — 


zd H 


-Xlo- 


Xle4 


Ylo- 


yie4 


ZI0 + 


Zle 


-X2 


-Y2 + 


Z2 


H4 = 


+za + 


zb- 


zc — 


zd H 


-Xlo- 


A:ie4 


YI0 + 


yie4 


Zlo- 


Zle 


-X2 


+ Y2- 


Z2 


H5 = 


+za + 


zb + 


zc + 


zd H 


'XI0 + 


Ariel 


YI0 + 


yie4 


ZI0 + 


Zle 


+ X2 


+ Y2 + 


Z2 


He = 


+za — 


zb- 


zc + 


zd -\ 


-XI0 + 


Ale- 


YI0 + 


yie- 


ZI0 + 


Zle 


+ X2 


-Y2- 


Z2 
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Hj = -~za -zh-zc-zd- Xlo + Xle - Ylo + Yle + Zlo + Zle - X2-Y2 + Z2 
Hs = -za + zb + zc- zd- Xlo + Xle + Ylo + Yle - Zlo + Zle -X2 + Y2- Z2. 

I 

The sum of these eight intervals lead to the cancellation teraction, so that H = Hnn/2 + 0{At). 
of the one-body terms and the next-nearest neighbor in- 
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